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ABSTRACT The DNA dannage response (DDR) is a dynannic process that is crucial for protecting the cell KEYWORDS 
from challenges to genonne integrity. Although many genome-wide studies in Saccharomyces cerevisiae genetic 
have identified genes that contribute to resistance to DNA-damaging agents, more work is needed to interaction 
elucidate the changes in genetic interaction networks in response to DNA lesions. Here we used conditional profiles 
epistatic miniarray profiling to analyze the genetic interaction networks of the DDR genes RTT107, SLX4, DNA damage 
and HRQ1 under three DNA-damaging conditions: camptothecin, hydroxyurea, and methyl methanesulfo- response 
nate. Rtt107 and its interaction partner Slx4 are targets of the checkpoint kinase Med, which is central to helicase 
the DDR-signaling cascades. Hrq1 recently was identified as a novel member of the RecQ helicase family in conditional 
S. cerevisiae but is still poorly characterized. The conditional genetic networks that we generated revealed interactions 
functional insights into all three genes and showed that there were varied responses to different DNA 
damaging agents. We observed that RTT107 had more genetic interactions under camptothecin conditions 
than SLX4 or HRQ1, suggesting that Rtt107 has an important role in response to this type of DNA lesion. 
Although RTT107 and SLX4 function together, they also had many distinct genetic interactions. In particular, 
RTT107 and SLX4 showed contrasting genetic interactions for a few genes, which we validated with indepen- 
dently constructed strains. Interestingly, HRQ1 had a genetic interaction profile that correlated with that of 5LX4 
and both were enriched for very similar gene ontology terms, suggesting that they function together in the DDR. 



Mapping of genetic interactions has been a valuable and powerful 
approach to reveal connections within complex biological systems 
(Baryshnikova et al 2013). Much of this work has been done in 
Saccharomyces cerevisiae because of the tools in place to create double 
mutants and the availability of vast arrays of mutant libraries. Al- 
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though providing great biological insights, most screens to date have 
been conducted under unperturbed growth conditions, whereas many 
networks in cells respond to environmental stimuli. 

A significant type of environmental stimuli is DNA damage, which 
can be caused by external factors, such as exposure to genotoxins or 
ultraviolet light, or internal factors, such as replication fork stalling or 
DNA polymerase error (Lindahl 1993). Cells are constantly exposed to 
these insults that, if not properly repaired, may compromise genomic 
integrity or ultimately lead to cell death. Because of the vital impor- 
tance of genomic integrity, cells have complex mechanisms to regulate 
the DNA damage response. DNA damage is detected by sensors, 
which trigger a signaling cascade, leading to the activation of the 
kinases Mecl and Tell, the yeast homologs of mammalian ATR 
(ATM and Rad 3-related) and ATM (ataxia-telangiectasia mutated). 
These kinases in turn elicit various cellular responses, including cell- 
cycle arrest, DNA repair, apoptosis, and/or DNA damage-induced 
transcriptional program (Putnam et al 2009; Finn et al 2012). 
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One of the downstream phosphorylation targets of Mecl is Rttl07/ 
Esc4, which is required for reinitiating rephcation after repair of alky- 
lating DNA damage (Rouse 2004; Roberts et al 2006). Deletion of the 
RTT107 gene results in hypersensitivity to DNA-damaging agents 
such as the DNA-alkylating agent methyl methane-sulfonate (MMS), 
the nucleotide reductase inhibitor hydroxyurea (HU), and the topoiso- 
merase I poison camptothecin (CPT) (Chang et al 2002; Rouse 2004; 
Parsons et al 2006; Roberts et al 2006). Rttl07 contains several BRCT 
(BRCAl C-terminal) homology domains, which often serve as phos- 
phobinding modules to recruit signaling complexes and repair factors to 
DNA damage-induced lesions (Rouse 2004; Mohammad and Yaffe 
2009). Consistent with a role as a scaffold for protein-protein interac- 
tions during the DNA damage response, Rttl07 interacts with a number 
of DNA repair and recombination proteins and is recruited to sites of 
DNA lesions (Chin et al 2006; Roberts et al 2006, 2008; Ohouo et al 
2010; Leung et al 2011; Ullal et al 2011). 

The best-characterized Rttl07-interacting partner is the replication- 
specific endonuclease Slx4, which interacts with the N-terminal BRCT 
domains of Rttl07 (Roberts et al 2006). Slx4 is required for Mecl- 
dependent phosphorylation of Rttl07 and, like Rttl07, facilitates resump- 
tion of DNA replication after DNA damage (Roberts et al 2006). However, 
it has become clear over the last few years that Rttl07 also has Slx4- 
independent functions, and vice versa. Consistent with this, the defects 
in DNA damage response are generally more severe in rttl07A mutants 
than in slx4A mutants, and rttl07A sbc4A double mutants are more sensi- 
tive to MMS than either of the single mutants (Roberts et al 2006). 

Although Slx4 has been studied in the context of its interaction 
with Rttl07, the SLX4 gene was first identified in a synthetic lethal 
screen with SGSl, which encodes for a RecQ helicase (Mullen et al 
2001). DNA helicases represent an important class of enzymes in- 
volved in the DNA damage response and have roles in recognition 
of DNA damage, DNA recombination, and stabilization of stalled 
replication forks (Brosh 2013). In S. cerevisiae, Sgsl was thought to 
be the only RecQ helicase family member until recently, when Hrql 
was identified as a novel RecQ helicase (Barea et al 2008; Kwon et al 

2012) . The functions of Hrql have only been preliminarily character- 
ized, but based on the relationship between Slx4 and Sgsl, Hrql may 
also have linkages to Slx4 that have yet to be uncovered. 

Initial genome-wide studies to characterize genetic function in 
response to DNA damage measured the fitness of deletion mutants 
exposed to a variety of genotoxic insults (Giaever et al 2004; Parsons 
et al 2006; Hillenmeyer et al 2008). However, these studies only 
evaluated the requirement of single genes for resistance to DNA- 
damaging agents, whereas the effects on genetic networks were only 
studied in a small directed screen (St Onge et al 2007). Two recent 
studies used genetic interaction mapping to gain new insights into the 
DNA damage response (Bandyopadhyay et al 2010; Guenole et al 

2013) . In the initial study, all possible double mutants of 418 genes 
were created and exposed to MMS to evaluate changes in the genetic 
interaction network (Bandyopadhyay et al 2010). Using this approach, 
the authors demonstrated that differential genetic interactions are bet- 
ter able to reveal functions in the DNA damage response and identified 
new roles for several genes. A follow-up study expanded on this work 
and interrogated 55 query genes crossed to a library of more than 2000 
genes in MMS, CPT, and zeocin conditions (Guenole et al 2013). 
Analysis of the differential genetic networks revealed several genes that 
were hubs of genetic interactions, and additional experiments demon- 
strated that these genes had novel roles in the DNA damage response. 

Here we use a similar approach of measuring conditional genetic 
interactions to study further the flmctions of Rttl07, Slx4, and Hrql. 
We analyzed the significantly interacting gene pairs to identify those 



that emerged or changed in response to DNA damage. Overall, 
RTT107 exhibited more genetic interactions than SLX4 or HRQl in 
CPT conditions, indicating an important role for Rttl07 in responding 
to CPT. Furthermore, SLX4 and RTT107 showed distinct, and some- 
times even opposing, genetic interactions, even though the protein 
products exist at least in part as a complex in the cell. Interestingly, 
the interaction profile and enriched gene ontology terms for HRQl 
most closely resembled that of SLX4, suggesting that they have over- 
lapping flmctions in the DNA damage response. 

MATERIAL AND METHODS 
Yeast strains 

AU yeast strains used in this study are listed in Supporting Information, 
Table SI and created using standard yeast genetic techniques (Ausubel 
1987). Complete gene deletions were achieved using one-step gene in- 
tegration of polymerase chain reaction-amplified modules (Longtine 
et al 1998). Mutants for conditional epistatic miniarray profiling (cE- 
MAP) screens were constructed in the BY4742 background, whereas all 
other strains were constructed in the W303-1A background. 

cE-MAP 

cE-MAP screens were performed and normalized as described 
previously (Collins et al 2010), with the exception that we expanded 
the number of plates per query strain to accommodate all drug con- 
ditions tested (see Figure 1). In brief, deletion mutants of the query 
genes were crossed, using a Singer robot, to a library of 1536 mutants 
covering a number of categories, including kinases/phosphatases and 
chromatin biology. We used 10 and 15 |jLg/mL of CPT, 50 and 
100 mM of HU, and 0.0075 and 0.0125% of MMS (Sigma- Aldrich), 
along with two no- drug controls. All strains and conditions were 
screened in triplicate. Complete cE-MAP profiles can be found in File SI. 

Differential S- scores were calculated by subtracting the S- score in 
yeast extract peptone dextrose from the S-score in each drug condition 
for each gene pair, and converted to Z-scores. Corresponding p-values 
were corrected for multiple testing using the fdrtool R package (Strimmer. 
2008) and conditional genetic interactions were considered significant for 
q-vdluQS < 0.05. cE-MAP networks were visualized with Cytoscape 
(Cline et al 2007). 

Growth and DNA damage sensitivity assays 

Overnight cultures grown in yeast extract peptone dextrose were diluted 
to 0.5 OD600. The cells were tenfold serially diluted and spotted onto 
solid yeast extract peptone dextrose plates or plates with MMS, CPT, or 
HU (Sigma- Aldrich) at the indicated concentrations. The plates were 
then incubated at 30° for 2 d and subsequently photographed. 

Gene Ontology (GO) analysis 

The Database for Annotation, Visualization, and Integrated Discovery 
was used for GO term enrichment analysis (Dennis et al 2003). For each 
query gene, a list of significantly interacting genes was compiled that 
included aU drug conditions and both positively and negatively interacting 
genes. Multiple testing correction was done using the Benjamini method, 
and enriched GO terms were considered significant for ^-values < 0.05. 

RESULTS 

Genetic interaction profiles were considerably altered 
when exposed to DNA-damaging agents 

To gain a global understanding of the functions of Rttl07, Slx4, and 
Hrql, we measured genetic interactions with representative genes 
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Double Mutant 



No Drug HU MMS CPT No Drug HU MMS CRT 

SOmM 0.0075% lOug/ml lOOnoM 0.0125«» 15ug/ml 

Figure 1 Schematic diagram of the conditional epistatic miniarray pro- 
filing (cE-MAP) workflow. The number of plates was expanded at various 
points in the process to accommodate for the multiple conditions tested. 
All screens were performed three times. HU, hydroxyurea; MMS, methyl 
methane-sulfonate; CPT, camptothecin. 

across the genome under three DNA damaging conditions (MMS, 
HU, and CPT). To achieve this, we used a version of the synthetic 
genetic array technology, the epistatic miniarray profile (E-MAP), to 
map genetic interactions for RTT107y SLX4, and HRQl under differ- 
ent DNA-damaging conditions, which we termed cE-MAP (Tong 
et al 2001; Baryshnikova et al 2010; Collins et al 2010). During 
the mutant selection process, the plates were expanded to accommo- 
date the four different conditions tested and to measure each condi- 
tion with three technical replicates (see Figure 1 for workflow). 

The analysis pipeline for E-MAP calculates S- scores, which reflect 
both the strength of the genetic interaction and the statistical con- 
fidence (Collins et al. 2010). For a broad assessment of how DNA 
damage affected the genetic networks, we first calculated the Pearson's 
correlations of the query genes' S-score profiles under the three drug 
conditions (Figure 2A). Supporting the idea that genetic networks 
respond significantiy to external stimuli, the genetic interaction pro- 
files generated under unperturbed growth conditions clustered away 
from the profiles generated under DNA-damaging conditions. Strik- 
ingly, the profiles for RTT107 in DNA-damaging conditions were 
closely correlated, regardless of the type of DNA insult. In contrast, 
the profiles for SLX4 and HRQl clustered together by the type of 
DNA-damaging agent. In addition to these major patterns, the genetic 
interaction profiles of all query genes under the same type of DNA- 
damaging agent also were positively correlated to one another (for 
example, see the positive correlations between the RTT107, SLX4, and 
HRQl profiles in MMS). Finally, we observed that the two concen- 
trations of drug used produced very similar profiles; thus, we averaged 
the scores of the two concentrations. This improved the confidence of 
our S-scores since they were represented by six replicates instead of 
three; these sets of average S-scores were used for all further analyses. 



According to previously published thresholds, S-scores lower than 
-2.4 or greater than 2.0 are considered significant genetic interactions 
(Collins et al 2010). Consistent with the dynamic nature of the DNA 
damage response, the distribution of the S-scores for each query gene 
changed significantly between the unperturbed growth conditions and 
the DNA-damaging conditions (Figure 2B and Figure SI). For exam- 
ple, under HU conditions, the S-score distribution for each query gene 
became broader, indicating that the number and strength of the ge- 
netic interactions increased (Kolmogorov-Smirnov test, two-sided, 
p- value < 1.5 X 10~^ for all comparisons). Strikingly, rttl07A mutants 
had the broadest distribution of S-scores of the three query mutants, 
i.e.y the most significant genetic interactions, followed by slx4A 
mutants, then hrqlA mutants (Figure 2B). We note that this was 
consistent with the relative DNA damage sensitivity of these mutants 
(Figure 2C). The distribution of the S-scores of all conditions for aU 
three query genes were significantiy different from one another, al- 
though the distributions for SLX4 and HRQl looked more similar to 
one another than to that of RTT107 (Kolmogorov-Smirnov test, two- 
sided, p-value < 7.8 X 10~^^ for all comparisons. Figure SID). 



RTT1 07 exhibited more conditional genetic interactions 
than SLX4 or HRQ1 

To identify genetic interactions under DNA damage conditions that 
were significantly different from unperturbed growth conditions, we 
adapted a published method (Bandy opadhyay et al 2010). Specifically, 
the differences between the S-score in the DNA damage condition and 
the S-score in unperturbed growth conditions for each gene pair were 
subtracted from the average of all the differential scores and divided 
by the standard error to calculate a Z score. Genes with significant 
Z scores after correcting for multiple testing are listed in Table 1 
(^-value < 0.05). In total there were 569 gene pairs found to have 
significant differential interaction in the drug conditions tested. Of these, 
378 were negative interactions (DNA damage-induced sickness or lethal- 
ity), and 191 were positive interactions (DNA damage-induced epistasis 
or suppression). To test the reliability of this approach, we looked for 
known condition-specific genetic interactions. Consistent with our pre- 
viously published results, deletion of DOTl and BREl suppressed the 
DNA damage sensitivity of rttl07A in MMS (Levesque et al 2010). 
Furthermore, the suppression by dotlA was additionally observed in 
HU but not in CPT, whereas the suppression by brelA was limited to 
MMS, in both the cE-MAP data and the independentiy constructed 
mutants (Figure 3). 

We visualized the conditional genetic interactions in a network in 
which nodes represented query or array genes and edges represented 
significant conditional genetic interactions. The edges were colored 
according to the drug condition that the genetic interaction occurred 
in (Figure 4A). As revealed by the network map, there was a subset of 
genes that interacted with all three query genes, suggesting that they 
play a more general role in the DNA damage response. These included 
the homologous recombination genes RAD52, RAD55, and RAD57 
(Figure 4B). Aside from this group of genes, there were also subsets 
that interacted with only two out of the three query genes. HRQl and 
SLX4 shared the greatest number of interacting genes, and this rep- 
resented a significant overlap between these two groups (Fisher's exact 
test, greater, p-value = 2.2 x 10~^^, Figure 4A). Further supporting 
shared functions of RTT107 and 81X4, there was also a significant 
overlap between their interacting genes (Fisher's exact test, greater, 
p- value = 1.7 X 10~^). Conversely, each query gene had unique genetic 
interactions, and RTT107 had the greatest number of these (Figure 
4A). Whereas the majority of the unique genetic interactions with 
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Figure 2 Genetic interaction profiles changed in response to DNA-damaging agents. (A) Pearson's correlation revealed that the genetic interaction profiles 
of RTT107were more similar to one another regardless of the drug condition, whereas SLX4 and HRQ1 were more similar to one another. (B) Density plots 
of the S-scores for each query gene showed a broader distribution of S-scores for RTT107 than SLX4 or HRQ1. (C) rtt107A mutants were more sensitive to 
DNA-damaging agents than s/x4A or hrq1A. 10-fold serial dilutions of the indicated strains were plated onto media containing the indicated drugs. 



RTT107 occurred under CPT conditions, HRQl and SLX4 had min- 
imal numbers of unique interactions in CPT. 

When comparing the total number of genetic interactions, we 
found that RTT107 had many more positively interacting genes than 
either SLX4 or HRQl (Figure 4C). After these interactions were sep- 
arated into each drug condition, it became clear that the biggest 
difference in interactions occurred during exposure to CPT for both 
positive and negative interactions (Figure 4, D and E). Taken together, 
these data suggest that Rttl07 plays an important role in responding 
to protein -bound nicks induced by CPT. To support this, RTT107 also 



had a strong positive/epistatic genetic interaction with TOPI, the 
molecular target of CPT, whereas this interaction was absent for 
SLX4 and HRQl (Figure 5A). 

Genes with condition-specific interactions were 
enriched for functions in the DNA damage response 

To further analyze the functions of RTT107, 81X4, and HRQl revealed 
by the conditional genetic interactions, we looked at the enrichment of 
GO terms using the Database for Annotation, Visualization, and 
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P Table 1 Genes that showed condition-specific interactions with query genes that were statistically significant after multiple test 
correction (q < 0.05) 



Query 



Drug 



Interaction 



Significant Genes 



HRQ1 CPT Negative ASF1, CHU, CLB5, CTF4, DCCI, DDC1, MMS1, MMS22, MRE11, PBY1, RAD17, RAD24, 

RAD52, RAD54, RAD55, RAD57, RAD59, RTT101, RTT109, SAE2 
HRQ1 CPT Positive CYC8, GMH1, PDA1 

HRQ1 HU Negative ARP4, ASF1, BAS1, BMH1, BRE1, CLB5, ERG5, ERJ5, GCN1, GCN20, GET2, GNP1, 

HPC2, IRA2, LAT1, LGE1, LST4, MET18, MFT1, MKS1, MRC1, MRE11, NPR1, PBS2, 
PDB1, PDE2, PMR1, POL32, RAD52, RAD54, RAD55, RAD57, RIM21, RPL34B, RPS21B, 
RTF1, RTG3, RTTW9, SDC1, SEC22, SEC66, SEC72, SGF73, SNX4, SPF1, SUA7, SWD1, 
SWD3, SWI4, SYC1, UBA3, UBP15, UBP6, URE2, VPS8, YTA7 

HRQ1 HU Positive AIM21, APS3, AR01, ARP7, BTS1, BUU, CAP1, COQ2, CUE3, CYC8, CYT1, DCR2, FEN1, 

HOS2, MSS18, NGL2, PET130, THIS, YMR102C 

HRQ1 MMS Negative DAL81, FKS1, GFD1, GFD2, GNP1, MMS2, MOG1, MPH1, MRE11, MSH4, PBY1, PET18, 

PH05, POL32, RADIO, RAD18, RAD27 , RAD59, REV1, REV3, REV7, RVS161, SAK1, SCS7, 
SR09, STP1, TRS33, UBC13, VMA21, YGL081W, YSY6 

HRQ1 MMS Positive CLB5, HST3, ILM1, KAP122, PMR1, SUR4 

RTTW7 CPT Negative AGE2, BCK1, BFA1, BUB2, CHS5, CLB2, CLB5, COG5, CRN1, CSG2, DCC1, DDC1, DEP1, 

ECM33, ELA1, ERVU, FET3, FKS1, FUN30, GAS1 , HPC2, IRA2, IRC21, KEX2, LAS21, LEM3, 
LGE1, LTE1, MAO 7, NCL1, OP/3, OST3, PAC1, PBY1, PEP8, PFA4, PMT1, PMT2, PPH21, 
PPH3, PPM 7, PRE9, PSP2, RAD17, RAD24, RAD54, RAD55, RAD57, RAD61, RDH54, RDI1, 
REV7, RGA1, RRD2, RTF1, RTS1, RXT2, SAP30, SCS7, SEC22, SEC28. SLA1, SM/7, SMy7, 
SPF7, SR09, STV1, SUR4, SWI4, TPK3, UBPU, VPS24, VPS27, VPS29, VPS35, VPS5, VPS8, 
VPS9, YDR061W, YJR088C, YLR426W, YPL150W, ZDS1 

RTTW7 CPT Positive APQ12, ARP7, BUB1, BUB3, CDC28, CYC8, CYT1, DBF2, DPB4, ELC1, ELG1, ESS1, GAC1, 

GAL80, HDA1, IKI3, IPT1, JHD2, LIA1, M/P7, MKK2, MMS4, MRC1, MRT4, NAP1, NCS2, 
NF/7, NOP12, OAF1, PAP2, PEF1, POL32, PSY4, RAD27, RPU1B, RPL8B, RPS1B, RPS4A, 
RSC4, SAC3, SAM37, SCT1, SFL1, SGF29, SLC1, SMT3, SPT2, SSN2, TED1, TOP1, TRM10, 
TUB3, TUP1, UBX4, UF01, YCR050C, YLR287C, YNR004W 

RTn07 HU Negative AGE2, BAS1, BUB2, CRN1, DUG2, GCN1, IRA2, MAD2, NPR1, OXR1, PBS2, PSP2, SEC22, 

SER2, SLX9, SLY41, SPF1, SUR4, SWI4, URE2, YDL089W 

RTT107 HU Positive AEP2, AIM26, AR01, ARP7, BTS1, BUB3, BUD21, BUD6, BUU, CDC48, CYC8, CYK3, CYT1, 

DBP1, ECM5, ERG6, HDA1, JNM1, LAG1, MET22, NAP1, NUP2, PEF1, PET130, PEX17, 
PMR1, POL30, PPM1, RIM101, RPN13, RPS1B, RPS28B, RPS30A, SCT1, SEF1, SPE3, SPT21, 
SSE1, THIS, TSA1, UBX4, UF01 

RTT107 MMS Negative BCK1, BFA1, CLB2, CRN1, DDC1, DUG2, FKS1, GFD1, GNP1, HRT1, IRC21, LSM7, LTE1, 

MBP1, MMS2, PAC1, POL30, PPM 7, PTC2, RAD17, RAD 18, RAD27, REV7, RTS1, SCS7, 
SLX9, SR09, SRS2, STP1, TEL1 , TSA1, UBC13, VMA21, YGL081W, YPL041C 

RTTW7 MMS Positive AIM29, APQ12, AR01, BMH1, BTS1, CSM3, CYT1, DBF2, DOT1, ELC1, FEN1, GAL80, GSF2, 

HST3, ILM1, LGE1, MRC1, MSS18, NAM7, NPP1, PER1, RAD52, RMD11, ROT2, RPL8B, RPS1B, 
RPS4A, RSC4, SPT2, TEP1, TMA23, TOP1, TUB3, UF01 

SLX4 CPT Negative ARP4, ASF1, CHU, CIK1, CLB5, CSM3, CTF4, DCC1 , DDC1, DIA2, GAS1, LAS21, LEM3, LGE1, 

MMS7, MMS22, PBy7, PMR7, PPH3, RAD17, RAD24, RAD52, RAD54, RAD55, RAD57, RAD59, 
RDH, RTT101, RTn09, SAE2, SGF73, SRS2, STV1 

SLX4 HU Negative AGE2, AIM32, ARL3, ARP4, BAS1, BMH2, BRE1, CLB5, COG5, CWH41, DDC1, ECM30, ERG5, 

ERJ5, FKH2, GCN1, GCN20, GEF1, GNP1, HCM1, HPC2, INP53, IRA2, LGE1, LST4, MDS3, 
MRE11, NPR1, OXR1, PBS2, PBY1, PDB1, PDE2, PFA4, PMR1, PRE9, PSP2, QCRW, QNQ1, 
RAD54, RAD55, RAD57, RPL41B, RPS11A, RTF1 , SDC1, SEC22, SEC66, SGF73, SKY1, SNG1, 
SPF1, STV1, SUA7, SUR4, SWD1, SWD3, SWI4, SYC1, SYF2, UBP15, URE2, VPS27, YDROSIW, 
YER064C, YPR063C, YSY6 

SLX4 HU Positive APS3, BTS1, BUU, CAP1, CDC36, COQ2, CYT1, FEN1, HOS2, IMP2, MDM38, MFT1, MSS18, 

NGL2, PET123, PET130, RPL43A, RPN13, SCD6, SPT21, TAF9, TOP1, YDU76W, YIP3 

SLX4 MMS Negative CKB1, DIA2, ENT4, ERV25, GFD1, GNP1, HST3, IMP2, MMS2, MPH1, PAC1, PET18, POL32, 

PPH3, PSY3, RAD18, RAD26, RAD27, RAD59, REV3, REV7, RRD1, SAE2, SAP185, SCS7, 
STP1, TOM7, TSR3, UBC11, UBC13, VMA21, YGL08m 

SLX4 MMS Positive BUD14, DOT1, MSH4, NAM7, RPS21B 



CPT, camptothecin; HU, hydroxyurea; MMS, methyl methane-sulfonate. 



Integrated Discovery (Dennis et al. 2003). For each query gene, all the 
genes that had significant conditional interactions in all three DNA- 
damaging conditions were analyzed for GO term enrichment. Because 
the library of genes tested with our cE-MAP approach was already 
enriched for nuclear function, we used this list of genes as the back- 
ground for the analysis rather than the whole set of genes in the yeast 



genome. AU the significantly enriched GO terms were related to the 
DNA damage response as expected (Table 2), indicating that the 
treatments induced DNA damage specifically and not a general stress 
response. However, the specific enriched processes differed between 
the three query genes. Both SLX4 and HRQl were enriched for pro- 
cesses related to DNA metabolism and repair, and postreplication 
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Figure 3 cE-MAP data recapitulated the drug-specific genetic inter- 
actions between RTT107 and DOT1 or BRE1. (A) Shown are subsets of 
cE-MAP data. Blue and yellow represent negative and positive genetic 
interactions, respectively. (B) 10-fold serial dilutions of the indicated 
strains were plated onto media containing the indicated drugs. 

repair was one of the most overrepresented terms (5.45- and 6.42-fold 
enrichment, respectively). Notably, the lists of enriched GO terms for 
SLX4 and HRQl were almost identical, further suggesting that these 
two genes have similar functions in the DNA damage response. In 
contrast, RTT 107 wdis significantly enriched only for one GO term, cell 
cycle checkpoint (2.95 fold enrichment, q < 0.03), although many of 
the GO terms that did not meet the significance cutoff were also 
related to the cell cycle (data not shown). 

Deletion of HST3 and MRC1 suppressed the DNA 
damage sensitivity of rtt107A but not slx4A mutants 

Hierarchical clustering of the significantly interacting genes revealed 
several patterns of genetic interaction profiles. There were sets of genes 
that were specific for the query gene, regardless of the DNA-damaging 
agent (Figure 5A). Intriguingly, there were certain genes that showed 
strong positive interactions with RTT107 but negative interactions 
with SLX4, such as HST3 and MRCl, further supporting the idea 
that Rttl07 and Slx4 have unique functions. 

Conversely, other sets of genes were specific for the DNA- 
damaging agent and interacted with all three query genes under that 
condition (Figure 5B). For example, the CPT-specific genes included 
RAD24, RAD17y and DDCU which encode for components of the 9-1- 
1 checkpoint clamp and RFC loader complex (Majka and Burgers. 
2007). The MMS-specific genes included multiple components of the 
translesion synthesis pathway, such as REV3, REV7, and RAD 18 
(Sharma et al 2013). Unexpectedly, the HU-specific genes included 
several transcription-related genes such as SWDl, SWD3, and SDCl. 

Positive S-scores indicate that the double mutant exhibits better 
fitness than expected (multiplicative product of single mutants' fitness), 
but it does not differentiate between suppression and epistasis. To 
further investigate the nature of the genetic interaction between 
HST3 and RTT107 and SLX4, we independently constructed deletion 
mutants and extended the analysis to include HST4, which was not on 
the E-MAP library. Hst3 and Hst4 are protein deacetylases that are 
both responsible for removing histone H3 K56 acetylation, thereby 
affecting replicative lifespan and response to DNA damage (Miller 
et al 2006). Deletion of HST3 clearly suppressed the DNA damage 
sensitivity of rttl07A mutants in all three drugs tested, albeit to a lesser 
extent in HU (Figure 6A). Confirming the striking opposite interac- 
tions of RTT107 and SLX4 observed in the E-MAP, HST3 and SLX4 



showed a synergistic interaction in CPT and MMS, but not HU. In 
general, HST4 showed the same genetic interaction profile as HST3. 
However, the hst3Ahst4Arttl07A triple mutant showed a variable phe- 
notype depending on the DNA-damaging agent, and this differed from 
the hst3Arttl07A or hst4Arttl07A double mutants, portraying a com- 
plex relationship between HST3 and HST4 in this genetic interaction. 

As a second example of characterizing positive interactions, we 
focused on MRCl, a gene encoding for an S-phase checkpoint adaptor 
(Tanaka 2010), which also showed this pattern of opposite genetic 
interactions with RTT107 and SLX4. Using independently constructed 
deletion mutants, we observed that deletion of MRCl suppressed the 
DNA damage sensitivity of rttl07A mutants to CPT and MMS (Figure 
6B). Interestingly, the converse was observed in the case of HU, in that 
deletion of RTT107 mildly suppressed the sensitivity of mrclA 
mutants. In contrast, the slx4AmrclA double mutant was clearly more 
sensitive to CPT and MMS than the slx4A single mutant, although the 
deletion of SLX4 mildly suppressed the HU sensitivity of the mrclA 
single mutant. 

DISCUSSION 

In this study we used cE-MAP to generate conditional genetic 
interaction profiles for RTT107y SLX4, and HRQl, to further investi- 
gate their functions in the DNA damage response. We tested CPT, 
HU, and MMS, which elicited specific genetic interactions with each 
of the query genes, and the network of interactions observed provides 
insight into the mechanisms of both the DNA-damaging agent and 
the query genes. Furthermore, we validated two specific examples of 
genetic interactions emerging from the cE-MAP in direct genetic tests 
and identified a novel genetic suppression in each case. 

A critical component of the cE-MAP approach was the use of 
DNA-damaging agents. Because we are interested in interrogating 
the DNA damage response functions of the query genes, it was 
crucial to evaluate the genetic interactions in conditions when 
those functions are active. Testing three different DNA-damaging 
agents also provided an opportunity to compare between the 
responses of the genetic interaction network to each type of DNA 
insult. 

We observed sets of genes that showed significant interactions 
under each specific drug condition. For each drug condition there was 
different DNA damage response pathways represented within the sets 
of genes, supporting the idea that the cell responds specifically to 
different types of DNA lesions. The CPT-specific genes included 
components of the 9-1-1 checkpoint clamp and RFC loader complex, 
suggesting that the DNA damage response to protein adducts involves 
this component of the checkpoint response (Majka and Burgers. 
2007). The MMS-specific genes included multiple components of 
the translesion synthesis pathway, which is one of the pathways in 
postreplication repair that allows cells to replicate past damaged bases 
or bulky adducts (Sharma et al. 2013). These data suggest that it also 
plays a role in bypassing alkylated bases. The HU-specific genes in- 
cluded several transcription-related genes, which was unexpected, but 
could be explained by an indirect effect of the mechanism of HU, 
which depletes the deoxynucleotide triphosphate pool. Interestingly, 
one of the HU-specific genes, YER064C, is relatively uncharacterized 
but was recently shown to change its cellular localization upon expo- 
sure to HU, suggesting a role for this gene in response to replication 
stress (Tkach et al 2012). Consistent with our study, the previously 
pubUshed E-MAP analysis of the DNA damage response revealed that 
genes showing significant interactions in CPT are enriched for func- 
tion in the DNA damage checkpoint, whereas significant genes in 
MMS are enriched for post- replication repair (Guenole et al 2013). 
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Figure 4 RTT107 had more significant genetic interactions than SLX4 or HRQ1. (A) Visualization of all the significant genetic interactions. Nodes 
represent query or array genes, and edges are colored by the drug condition that the interaction occurred in. Blue represents CPT, green 
represents HU, and red represents MMS. The numbers of unique interactions for each query gene are labeled for each drug condition. Indicated 
p-values are from Fisher's exact tests (greater) of the genes that interact only with two out of the three query genes. (B) Enlarged view of a subset 
of the network indicated by the black box in (A). (C) RTTW7 had more positive genetic interactions than SLX4 or HRQ1. RTT107 had more (D) 
positive and (E) negative genetic interactions in CPT than SLX4 or HRQ1. CPT, camptothecin; HU, hydroxyurea; MMS, methyl methane-sulfonate. 
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Figure 5 Different patterns of genetic interactions were observed for genes that significantly changed their interactions in response to DNA- 
damaging conditions. Shown are subsets of conditional epistatic miniarray profiling (cE-MAP) data. Blue and yellow represent negative and 
positive genetic interactions, respectively. (A) Some genetic interactions were specific for the query genes. (B) Other genetic interactions were 
specific to the drug condition and were common across all query genes. 



In contrast to the drug- specific genes, there were also sets of genes 
that showed unique interactions with each specific query gene. These 
interactions suggested that the query genes we interrogated have 
distinct fianctions in the DNA damage response. It is of particular 
interest that RTT107 and SLX4 shared only a subset of genetic inter- 
actions, given that the Rttl07 and 81x4 proteins exist as a complex in 
the cell (Roberts et al 2006). An attractive model is that there are 
different pools of Rttl07 and Slx4 protein complexes that contribute 
to specific functions, since Rttl07's interaction with SMC5/6 and 
Slx4's interaction with Ski are independent of each other (Roberts 
et al 2006; Leung et al 2011). The human homologs of Rttl07 and 
Slx4, PTIP and SLX4, respectively, also have many distinct functions. 
While PTIP is involved in the DNA damage signaling cascades and 
DNA repair pathway choice (Gong et al 2009; Wu et al 2009; Callen 
et al 2013), SLX4 has roles in HoUiday junction resolution and telo- 
mere length regulation (Castor et al 2013; Garner et al 2013; Wan 
et al 2013; Wilson et al 2013; Wyatt et al 2013). The data from this 
cE-MAP provide an opportunity to further elucidate the unique func- 
tions of Rttl07 and Slx4, which may be further dissected into 
responses to different DNA lesions. 



The cE-MAP data also provided more insight into the fianction of 
the helicase Hrql, which has only been preliminarily characterized 
(Kwon et al 2012; Choi et al 2013; Bochman et al 2014). Interestingly, 
the genetic interaction profile of HRQl correlated more closely to SLX4 
than RTT107y and the sets of genes that interacted with SLX4 or HRQl 
overlapped significantly. In addition, GO analysis of the significantly 
interacting genes returned almost identical GO terms for HRQl and 
SLX4. The close relationship between HRQl and SLX4 revealed by the 
cE-MAP data are supported by previous studies showing that SLX4 is 
synthetic lethal with SGSl, the major RecQ helicase in S. cerevisiae 
(Mullen et al 2001). Moreover, there is some evidence suggesting that 
Hrql and Slx4 are both involved in interstrand crosslink repair and 
suppression of telomere addition (Zhang et al 2006; Ward et al 2012; 
Bochman et al 2014). However, sk4A and hrqlA mutants display 
different sensitivities to DNA-damaging agents, indicating they fianction 
separately as well (Flott and Rouse. 2005; Choi et al 2013; Bochman 
et al 2014). Further experiments are needed to determine the roles of 
Hrql and Slx4 in DNA structure maintenance and the nature of their 
relationship in these fianctions. Possible routes of inquiry can be sug- 
gested by additional examination of the genetic data. 
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W Table 2 GO terms that were significantly enriched in the list of genetically interacting genes under all conditions 
using DAVID 



QuGry 


GO TGrm 


Fold Enrichment 


C]f-value (Benjamini) 


RTTW7 


Cell-cycle checkpoint 


2.95 


2.56E-02 


SLX4 


DNA metabolic process 


2.54 


6.22E-07 




Response to DNA damage stimulus 


2.53 


1.87E-05 




DNA recombination 


4.09 


1 .29E-05 




DNA repair 


2.59 


4.70E-05 




Recombinational repair 


4.91 


6.39E-04 




DNA replication 


3.03 


7.59E-04 




Postreplication repair 


5.45 


1 .94E-02 




Cellular response to stress 


1.78 


2.31 E-02 


HRQ1 


DNA metabolic process 


2.83 


8.29E-08 




DNA recombination 


4.64 


7.05E-06 




Response to DNA damage stimulus 


2.77 


5.32E-06 




DNA repair 


2.90 


8.83E-06 




DNA replication 


3.44 


2.28E-04 




Cellular response to stress 


2.08 


9.93E-04 




Postreplication repair 


6.52 


5.21 E-03 




Recombinational repair 


4.89 


5.12E-03 




DNA-dependent DNA replication 


3.91 


6.35E-03 




Double-strand break repair 


3.67 


1.13E-02 




Double-strand break repair via single-strand annealing 


7.34 


2.97E-02 



GO, Gene Ontology; DAVID, Database for Annotation, Visualization, and Integrated Discovery. 



We followed up on two genes that had particularly striking 
conditional genetic interactions. Both HST3 and MRCl showed a strong 
positive genetic interaction with i^TTi 07 but a negative interaction 
with SLX4. Although Rttl07 and Slx4 form a complex, these ge- 
netic interactions suggest that not only do Rttl07 and Slx4 have 
independent functions, but they may have opposing functions in 
these contexts. Using a direct genetic test, we found that deletion of 
HST3, as well as HST4, suppressed the DNA damage sensitivity of 
rttl07A mutants. The known target of the Hst3 and Hst4 deacety- 
lases is H3 K56 acetylation (H3 K56ac). Whereas deletion of HST3 
alone causes an increase in H3 K56ac, deletion of HST4 alone does 
not change the acetylation levels, and only in the double mutant 
are all H3 molecules completely acetylated (Celic et al. 2006). In- 
triguingly, the suppression of the rttl07A mutant phenotype in 
CPT was observed upon deletion of HST3 or both HST3 and 
HST4, but not HST4 alone. In contrast, deletion of either HST3 or 
HST4 alone was sufficient to suppress the DNA damage sensitivity of 
rttl07A mutants to MMS and HU. Based on this data, we speculate 
that the deacetylation of H3 K56ac may be important to the genetic 
interaction in CPT, whereas it is not relevant in MMS or HU con- 
ditions, rather there may be a different function or target of Hst3 and 
Hst4 involved. Similarly for the genetic interaction with 81X4, deletion 
of HST3 or HST4 alone exhibited the same phenotype, thus suggesting 
that deacetylation of H3 K56ac is not involved. 

We also validated the genetic interaction with MRCl and found 
that deletion of MRCl suppressed the DNA damage sensitivity of 
rttl07A mutants in CPT and MMS but aggravated the sensitivity of 
slx4A mutants. Interestingly, the situation was different in HU, where 
deletion of either RTT107 or SLX4 mildly suppressed the sensitivity of 
mrclA mutants. This finding is consistent with a model proposed by 
a previous study suggesting that Rttl07 and Slx4 inhibit the check- 
point adaptor protein Rad9, which is normally not important in rep- 
lication stress, but becomes crucial in the absence of Mrcl (Ohouo 
et al 2013). However, this model does not explain the genetic inter- 
actions observed in CPT and MMS and reflects the distinct responses 
to various types of DNA lesions, as well as the multiple functions of 
DNA damage response proteins. 



Our study contributes to the growing of body of data that has 
mapped genetic interactions in response to DNA damage and further 
validates it as a fruitful approach that reveals condition-specific functions 
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Figure 6 Positive S-scores from the cE-MAP data were based on 
suppression of DNA damage sensitivity of rtt107A mutants. 10-fold 
serial dilutions of the indicated strains were plated onto media con- 
taining the indicated drugs. (A) Deletion of HST3 and HST4 sup- 
pressed the DNA damage sensitivity of rttW7A mutants but 
aggravated the sensitivity of s/x4A mutants. (B) Similarly, deletion of 
MRC1 suppressed the DNA damage sensitivity of rtt107A mutants but 
aggravated the sensitivity of s/x4A mutants. 
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and pathways in the cell. There remains much ground to be covered as 
we have only started to characterize the pathways specific for the 
multitude of environmental conditions that affect all living organisms. 
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